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Abstract 

In the paper we study stochastic convolution appearing in Volterra equation 
driven by so called Levy process. By Levy process we mean a process with homo- 
geneous independent increments, continuous in probability and cadlag. 



1 Introduction 

Let if be a real separable Hilbert space with an inner product (•, and a norm | • |#. 
In the paper we consider a stochastic version of linear, scalar type Volterra equation in 
if of the form 

u(t)= / a{t-T)Au(T)dT + x + g{t), t > 0, (1) 
Jo 

where a G L 1 1 oc (lR + ), A is an unbounded linear operator in if with a dense domain T>(A), 
g is an if- valued mapping and x G H . 

The linear integral equation (JTJ is a subject of many papers connected with appli- 
cations in different fields. Among others, the equation may be applied to several 
problems arising in mathematical physics. For instance, theory of viscoelasticity provides 
numerous problems leading to the Volterra equation of the form (JTJ) (see [9], for survey). 
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We assume that the equation (JTJ) is well-posed and denote by TZ(t),t > 0, the family 
of resolvent operators corresponding to (JTJ) - Operators 1Z{t) are linear for each t > 0, uni- 
formly bounded on compact intervals, lZ(0)x = x holds on T>(A), and lZ(t)x is continuous 
on R + for each x G T>(A). Additionally, the following resolvent equation holds 



for all x G T>(A), t > 0. For more details concerning well-posedness and resolvent opera- 
tors we refer again to the monograph [9]. 

In the paper we study equation (0) with an external force g(t) = Z(t),t > 0, where 
Z is a Levy process defined on a stochastic basis (Q, J 7 , T t) P). This way the traditional 
Gaussian framework, when a Wiener process is the external noise, is extended. Our 
considerations are motivated by the growing interest in Levy processes in applications, 
when the empirical observations simply cannot be explained by means of the Gaussian 
distribution. 

So, we arrive at the equation 



Stochastic equations of Volterra type driven by semimartinagales have been of course 
already studied by some authors, for instance [7] or [11, 12]. But our paper treats the 
subject in a different spirit. We use the resolvent operators of the equation considered, 
then this way we try to extend the semigroup approach to the equation (j2J). For defining 
the stochastic convolution we do not use the general semimartingales technique but a 
simpler method. 

We have assumed that A is a closed linear operator in H with the dense domain T>(A). 
For such a class of operators exists the family TZ(t), t > 0, of resolvent operators for (j2J) 
which fulfills some useful properties (see again [9], Chapter 1). Moreover, the family 
lZ(t), t > 0, is a good enough class of operators to be integrands in stochastic integral 
with respect to a process with independent increments. 

Now, we introduce definitions of solutions to 0, analogously like in previously con- 
sidered stochastic cases. 

Definition 1 An H-valued predictable process X(t), t G [0, T], is said to be a weak 
solution to HJJ ; if P(j*\a(t - r)X(r)\ H dr < +oo) = 1 and if for all £ G V(A*) and all 
t > the following equation holds 





where X G H. 



(2) 



(X{t),0 H = (X , 0h+([ a(t - r)X(r)dr, A*£) H + (Z(t),£) 



H > 



P-a.s. 
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Definition 2 An H-valued predictable process X(t), t £ [0, T], is said to be a mild 
solution to if P(J |X(r)|^<ir < +oo) = 1 and, for arbitrary t £ [0, T\, 

X(t) =TZ{t)X + [ n(t-T)dZ{r), (3) 
Jo 

where 7Z(t) is the resolvent for the equation 

The aim of the paper is to study process called stochastic convolution 

Z n {t):= [ K{t-T)dZ(r), t>0, (4) 
Jo 

which is the crucial part of the solution ©. 

In section 2 we recall the rigorous definition of stochastic convolution with Levy pro- 
cess. Then we adapt it to the Volterra equation (J2J) driven by Levy process and use some 
properties of such integral. Next, in section 3 we consider particular Volterra equations. 

2 Stochastic convolution 

Assume that Z(t), t > 0, is an H- valued process with homogeneous independent incre- 
ments (that is, Levy process), defined on a fixed probability space (Q, J 7 , Tt-, P), contin- 
uous in probability, cadlag and Z = 0. It is of great importance in the study of linear 
and nonlinear stochastic Volterra equations driven by Levy processes to establish first the 
basic properties of the process Z-ji(t) = f* 1Z(t — r) dZ{r). 

The stochastic convolution (J3J), where Z is a Levy process, may be defined analogously 
like stochasic integral in the paper [2], that is, as a limit in probability of Stieltjes sums. 
This integral coincides with the integral defined by the semimartingales technique (see [6] 
or [8]), but approch used in [2] provides immediately some useful properties of the integral. 
The most important is that we obtain the explicit formula for the characteristic form of 
the convolution (jlj). 

First we recall some facts concerning stochastic integral with respect to Levy process 
used in the paper. 

The class C^ u , (H, G) of integrands is defined as follows: 
$ : [u, w] — > L(H,G), where u,w £ K + and H,G are Hilbert spaces, such that: 

1) for any h £ H, $/i : [u, w] — > G is measurable, 

2) \ \Q(s)\\ 2 ds < +oo, where || • || means the operator norm. 

(In the above definition L(H, G) denote the space of linear bounded operators acting from 
H into G.) 
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Theorem 1 (Theorem 1.3, [2]) 

Let H and G be Hilbert spaces and a function $ belong to the class C?JH, G). Assume 
that there exists a sequence {$„} of step functions that: 

• for any n G iV , $ n : [u, w] — ► L(H, G) and for partition u = Sq < s± < . . . < s n = t, 
$ n (s) = where s G (ss k+1 ], k = 0, 1, . . . , n - 1; 

• t/je following condition holds 

$(s)/i = lim $ n (s)/i /or any h E H for a. a. s E [it, w]; (5) 

n— >oo 

• t/iere exists a function g G L 1 ([n, w]) i/jai 

sup ||$ n (s)|| 2 < o(s) for a. a. sE[u,w]. (6) 

T/ien i/ie sequence of random variables 

„ n— 1 

J(* B ) := / S n (a) dZ s = £ *£(Z Sfc+1 - 

■/(w,to] fc=0 

converges in probability and the limit does not depend on the choice of the sequence {$ n }- 
Theorem 2 (Theorem 1.8, [2]) 

Let the function $ and the sequence {& n } satisfy the assumptions of Theorem^ Then 
the integral defined as 

[ $(s)dZ s :=P- lim f ® n {s)dZ s (7) 

J(u,w] n ^°° J(u,w] 

is well-defined G-valued random variable which has infinitely divisible distribution. 

Remark 1 The integral JJ° dZ s is defined as the limit in probability, as w — > oo ; 
of the integrals f,,Q(s)dZ s . The integrals f^ oo ^(s)dZ s and f*™ <&(s) dZ s are defined 
analogously. 

As we have already written, the aim of the paper is to study the stochastic convolution 
(jlj), where Z(t), t > 0, is a Levy process and TZ(t), t > 0, are the resolvent operators 
to the Volterra equation (j2J). Basing on properties of resolvent operators we can see that 
the operators TZ(t), t G [u, w], belong to the class C? W ,(H, H). Actually, TZ(t) are linear 
and bounded for each t > and lZ(t)x is continuous on M. + for each x belonging to the 
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domain T>(A) of the operator A. Moreover, the function !Z(-)x, x G T>(A), is measurable. 
Let TZ n (t), t G [u, w], be step functions defined as follows 

n 

ftn-Z^MXfc-!*], (8) 
i=l 

where t = u, to < t\ < . . . < £j_i < U < . . . < t n = w, and Sj is a point from [ii_i, ij]. 

Let us notice that the functions 7£ n defined by (JHJ) on the interval [m, iw], satisfy con- 
ditions (0) and (JBj). Indeed, the interval [u,w] is a compact set in K and the operator 
72.(t), t G [«, w] is continuous with respect to t. Then the sequence (JZ n ), n G AT, of step 
functions (JHJ is uniformly convergent to the function TZ. Additionally (7Z n ), n G N, are 
bounded. 

So, we may define the stochastic convolution (J1J) like the stochastic integral (|ZJ), that 
is, like the limit in probability of integrals of step functions. Hence, the following theorem 
comes directly from Theorem El 

Theorem 3 LetlZ(t), t > 0, be the family of resolvent operators of the Volterra equation 
(HJ) and Z(t), t > 0, be a Levy process. Then the integral 

[ K{t-T)dZ{r) := P - lim / TZ n (t-r)dZ(r) 
is well-defined H -valued random variable which has infinitely divisible distribution. 

Let us recall that the process Z(t), t > 0, as the process with independent increments, 
has the following representation (see for instance [1,3,5] or [10]) 

Z t = at + W t + A t , 

where a G H, (Wt),t > 0, is an if-valued Wiener process and (A t ),t > 0, is a jump process 
independent of (W t ). This decomposition is clearly unique. Moreover, for any t > 0, the 
random variable Z t has Levy characterization [ta,tO,tM], where is the covariance 
operator of W\ and M is the Levy spectral measure of Ai. This is a consequence of the 
fact that any infinitely divisible probability measure can be viewed as the distribution of a 
Levy process evaluated at time 1 and vice versa. Particularly, the famous Levy-Khintchine 
formula determines the class of characteristic functions corresponding to infinitely divisible 
laws. 

Now, we are ready to characterize the convolution (|3J) as follows. 
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Theorem 4 Let Levy process Z(t), t > 0, be such that every random variable Z t has 
Levy characterization [ta,tQ,tM]. Then the stochastic convolution 



Proof: By Theorem El the stochastic convolution given by the formula Zn(t) = J Q 7Z(t — 

r) dZ(r) , t > 0, has infinitely divisible distribution. In order to provide the Levy 
characterization of Zn(t) it is enough to write the characteristic functional of the law 
of Zfc(t) and next use the Levy-Khintchine formula for the corresponding characteristic 
exponent. The proof of the theorem is analogous to the proofs of Lemma 1.5 and Theorem 
1.8 from the paper [2]. In our case, the characteristic functional is 

f(y) := exp f Q (p(JZ(t — s)y)ds, for y G H, where (f)(w) = logE(expi(w, Zi)), w E H. □ 

Corollary 1 The stochastic convolution ^ is stochastically continuous and then has a 
predictable version. 

Theorem 5 Assume that the operators TZ(t), t > 0, are as above. Then the stochastic 
Volterra equation fijty has exactly one mild solution. 

Comment: Theorem 03 comes from uniqueness of the resolvent TZ(t), t > 0, for determin- 
istic Volterra equation ([!} and the existence of mild solution to (fT|). 

Theorem 6 Assume that the operators A,TZ(t) and the process Z(t), t > 0, are like 
above and the function a 6 W / 1 J ) ' C 1 (]R + ). Let X be an H-valued predictable process with 
integrable trajectories. If for any t G [0, T] and £ G T>(A*) the equality 




o Jo Jh 



Q 




(9) 



M 




(X(tU)H 



[ (a(t-T)X(r),A*C) H dT+ [ (^,dZ(r)) H 



(10) 



Jo Jo 



holds, then 



(11) 
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Proof: The idea of the proof is the following. First, we prove that if (110)1 is satisfied, 
then for any £ G C 1 ([0, T],V(A*)) and t G [0,T], the following equality holds 



*((d * X)(r) + a(0)X(r), rfr + /"(f(r), dZ(r)) H 



+ / (X(t),£(t))h, P-o.s. 



(12) 



Next, we take £(t) := - r)£ for r G [0,t] and rewrite JUJ). 

Then, using properties of resolvent operators, particularly the resolvent equation, we 
obtain thesis (JTTJ). □ 

Comment: The above Theorem |H1 says that a weak solution to (0) is a mild solution to 



Theorem 7 The stochastic convolution Z R (t) = J" 7£(t — r)dZ(r), t G [0,T] ; fulfills the 
equation HI]) . 

Proof: Let us notice that the process Zr has integrable trajectories (the set of ,, discon- 
tinuity" is at most countable). 
For any £ G T>(A*) we may write 

/ (a(t-T)Z R (T),A*Z) H dT = 
Jo 

(from®) = f (a(t-r) f K(r - a)dZ(a), A*£) H 
Jo Jo 



(from Dirichlet formula and stochastic Fubini theorem) 

rt r rt 



dZ(a),A*0n 



[ ( f a(t - t)K{t - a)dr 
Jo VJo 

rt V r-t—cr 

( / a(t - a - z)K(z)dz dZ(a),A*£) 
Jo Uo 



A[(a*K)(t-a)}dZ{a),OH 



(from resolvent equation) = ( / [lZ(t — a) — I]dZ '(a) , £) h 

Jo 



(/ K{t-a)dZ{a),0 H -{ dZ(a),0n 
Jo Jo 
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Hence, we obtained the following equation 

(Z R (t),£) B = ! (a(t-r)Z R (r),A*0 H dr+ [ (Z,dZ(r)) H 



o 



for any f 6 □ 
Corollary 2 Assume that the operator A is bounded. Then 



Z R (t) = [ a(t- T)AZ R (r)dr + [ dZ{r) . 
Jo Jo 



3 Particular cases 

Let us notice that till now we have not assumed that the resolvent operators TZ(t), t > 0, 
to the Volterra equation (J2J) have bounded variation. The stochastic integral with respect 
to Levy process and the stochastic convolution (jlj) have been defined and characterized 
(formula without this assumption. Hence, if the stochastic Volterra equation (J2J) 
driven by Levy process Z(t), t > 0, is well-posed, that is, the equation (J2J) has the 
resolvents TZ(t), t > 0, then Theorem |U for the convolution Z R (t), t > 0, holds. Besides, 
there are some Volterra equations which admit resolvent operators with bounded variation. 
These equations may be treated in a different way. So, in this section we study such 
particular case of Volterra equations. 

We will take the following common assumption. 

Assumption (A) There exists the family 7Z(t), t > 0, of resolvent operators to the 
Volterra equation (j2J), the operators have bounded variation and Z(t), t > 0, is an H- 
valued Levy process. 

If the resolvent operators TZ(t), t > 0, have bounded variation, we may use for stochas- 
tic convolution (J3J) the classical integration by parts. Then we may write 



TZ{t — r)dZ{r) — TZ(t — b)Z{b) — TZ{t — a)Z(a) — \ Z(r-)dK(t - r) , (13) 

(a,b] J(a,b] 

where Z(t), t > 0, is a stochastic Levy process. 

The properties of the stochastic integral J^ a fe , lZ(t — t)cIZ(t) may be obtained from the 
right hand side of (fT3*j) . Among others, we can deduce the following results. 

Proposition 1 Under the assumption (A) we have 

(x*, I n(t-r)dZ(r)) = [ n{t-r)d{x*,Z{r)) 

J(a,t] J(a,t] 

for x* G H* , where H* is the dual space to H . 
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Proposition 2 If the assumption (A) holds, the function 



t — ► / n{t - r)dZ{r) 

J(a,t] 

is H-valued random variable with infinite divisible distribution. 

This fact follows from (jlHjl and the approximation by Riemann-Stieltjes sums. See, 
e.g. [4]. 

Proposition 3 Let the assumption (A) be satisfied. Then 



log 



£( / K(t - r)dZ(r))(X)] = [ \og[C(Z(l))(\K(t - s))]ds , (14) 
J(o,t] J J(0,t] 



for A G H* , where H* is the dual space to H and C denotes the characteristic functional 
of the probability distribution of the appropriate random variable. 

Proof: The above formula (pPfj) comes from Lemma 1.1, [4] and the definition of the 
convolution (HI). □ 



Proposition 4 If the assumption (A) holds, then j^ ot ^IZ(t — r)dZ{j) converges 



in norm 



w.p. 1, as t — > +oc, if and only if L Qt ,TZ{t — r)dZ{r) converges in the distribution as 
t — > oo. 



Proof: Proposition comes from Lemma 1.2, [4]. □ 
We finish the paper by the examples of Volterra equations. 

Let us consider the case when the function a is completely positive. This class of 
kernels is very important in the theory of Volterra equations and arises naturally in 
applications. If the function a G L 1 1 oc (lR + ) is completely positive, then s(-,7), the solution 
to the integral equation 

s(t)+7 I a(t - t) s(t) dr = 1, t>0, (15) 
Jo 

is nonnegative and nonincreasing for any 7 > 0. More precisely, under this condition 
s(t) G [0, 1]. There is a relationship between the resolvent operator TZ(t) to the equation 
((21) and the corresponding function s(t, 7) fulfilling (|15)) . Namely, if —7 is an eigenvalue 
of A with eigenvector z 7^ 0, then lZ(t)z = s(t, j)z, t > 0. 
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Because every function monotonic on the interval [0,T] has bounded variation on 
[0,T], it is enough to choose any function a(t), t > 0, which is completely positive. 

Particularly, let H = L 2 (0, 1) and Au = D 2 u with V(A) = H 2 (0, 1) n ^(0, 1). The 
functions e&(£) = v/2/-7r sin k£, £ G [0,1], k G N, form an orthonormal sequence of 
eigenfunctions of the operator A, corresponding to the eigenvalues — /!& = —ir 2 k 2 , k G N . 
When we set additionally, that the function a(t) = e~* for t > 0, then we obtain in this 

case 

8(t,fi) = (1 +/ie- (1+,l)t ], t,fi > . 

There exists the resolvent TZ(t), t > 0, to the equation (J2J) in this case and is deter- 
mined by 1Z(t) Ck = s(t, Hkje-k, k E N . So, the function TZ(t), t > 0, is monotonic in 
[0, 1], because the function s{t,n) is. Hence, TZ(t) has bounded variation on [0, 1]. 

Other examples of Volterra equations with resolvent operators having bounded varia- 
tion may be found in the monograph [9] . 
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